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^ , Abstract 

1^-^ I We calculate the specific heat of the inner crust of neutron stars within a local-density approxima- 

f*"^ . tion to an improved relativistic Hartree-Bogoliubov theory. Non-uniformness of the system enhances 

' the specific heat in particular at low temperatures. The degree of enhancement is similar to that in 

, the spherical phase of Elgar0y et al. We examine a schematic interpolation between the results of 

Broglia et al. adopting the Gogny force and ours based on the Lagrangian of the relativistic mean 
' field model. 

> . 

The inner crust of neutron stars, with densities ranging from the neutron drip density to that of the 
^ , order of the saturation one, is a non-uniform system consisting of a Coulomb lattice of neutron-rich nuclei 

and a sea of free neutrons. This part conveys temperature drop due to neutrino emission in the core 
to the surface. The rate of the conveyance depends on the specific heat and the thickness of the inner 
crust. Since the former is sensitive to superfluidity, we study this non- uniform system as an application 
of the improved description of pairing correlation based on the RMF interaction, that is, a one-boson 
;h ' exchange interaction derived from the Lagrangian of the relativistic mean field (RMF) model. As a first 

step of this application, we adopt a local-density approximation and neglect possible non-spherical shapes 
of neutron-rich nuclei, in this talk. In this sense, this is a relativistic counterpart of Ref.Q and of the 
spherical part of Ref . . 

The first attempt to describe superfluidity in uniform infinite nuclear matter adopting the RMF 
interaction was done in Ref.l^; but its result was not very realistic in that the resulting pairing gap 
was too large. The reason can be ascribed to an unphysical behavior of the RMF interaction at high 
momenta: Repulsion is too strong. Note that repulsion also contributes positively to pairing gap at high 
momenta beyond a sign change of the gap function. Tanigawa and one of the present authors proposed a 
method to determine quantitatively a density-independent parameter to cut the momentum integration 
in the gap equation jj]. Although a more natural way to modulate the high-momentum part smoothly 
was also developed later [||, the simple method in Ref.Q is enough for the present purpose since only 
the pairing gap at the Fermi surface at each density is necessary under a local-density approximation. 

In order to describe the non-uniform system composed of a lattice of neutron-rich nuclei and a sea of 
free neutrons, we adopt the Wigner-Seitz (WS) approximation; each spherical cell contains at its center 
one nucleus permeated by free neutrons which dripped off from it. We parametrize the neutron density 
given by a density-dependent Hartree-Fock calculation in terms of a Woods-Saxon shape, 

Pn{r) = — r/^°'" p . , T + Poxt,« , (1) 

1 + cxp[(r - Rnj/an\ 
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as done in Ref. 0|. We perform calculations at five representative densities conforming to this reference; 
the parameters entering into Eq.(P and the radii of the WS cell are summarized in Table |l|. Here we 
introduce a local-density approximation — assume uniform Fermi gas of neutrons at each spatial point. 
Then the local Fermi momentum is calculated from the local density given by Eq.(pl) as 



fcF(r) = (37r2p„(r))i/3 



(2) 



Its profile in each WS cell is shown in the left panel of Fig.|i|. This determines the profile of pairing gap 
since the gap in the uniform neutron matter can be calculated as a function of by using the method 
of Ref. [Q. (Note that the symmetric nuclear matter case was presented there.) The result is shown in 
the right panel. The local pairing gap shown there gives the local quasiparticle energy, 



E{r,p) - {{Eip) - E{k^{r))f + A(r,p)2)i/2 ^ 
with E{p) = (p^ + M*^)-'^/^, M* being the density-dependent effective neutron mass. 







po,«[fm"^] 


Poxt,n[fm"^] 


i?n[fm] 


a„ [fm] 


i?ws[fm] 


(a) 


0.46 


0.114 


0.0737 


5.0 


1.2 


15.0 


(b) 


0.28 


0.101 


0.0436 


7.0 


1.1 


20.9 


(c) 


0.12 


0.098 


0.0184 


7.4 


0.8 


29.4 


(d) 


0.034 


0.104 


0.0047 


6.8 


0.9 


35.5 


(e) 


0.005 


0.108 


0.0005 


5.9 


0.9 


45.7 



(3) 



Table 1 : Parameters entering into Eq. (|^) and the radii of the WS cell. 
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Figure 1: Profile of the local Fermi momentum (left) and of the local pairing gap (right) in a Wigner-Seitz 
cell. Five curves in the order from the top to the bottom in the left panel correspond to densities (a) - 
(e) in Table [l[ 

Now we are ready to calculate the specific heat of the system. The specific heat is calculated from 
the entropy, 



S = -2 d 



_fp_ 
(27r)3 



{/(r,p) In f{r,p) + [1 - f{r,p)] ln[l - f{r,p)]} , 



Cvn - T— 



The phase-space distribution function in Eq.(|4|) is 



l + exp[S(r,p)/r] ' 



(4) 
(5) 

(6) 
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where the quasiparticle energy is given by Eq.(|). The final expression is 



Cvn — 



dr 



dp 



E{r,p) 



cos\i[E{r,p)/2T] 



E{r,p) 
T 



dE 

It 



(7) 



The spatial integration is performed over each WS cell and the momentum integration is done up to A = 
3.60 fm^^ determined in Ref.Q to give physical pairing gaps. Note that the temperature dependence of 
the quasiparticle energy is neglected because only low temperatures are of interest in the present study. 
In addition to neutrons, degenerate electron gas gives a comparable contribution, 



Cve = hlim^ + klY'^T , 
while that from protons is negligible. Then the total specific heat is given by 

Cytot = Cvn + Cve ■ 



(8) 



(9) 



The effect of non-uniformness can be seen by comparing the left (showing Cvn for uniform matter) and 
the center (for non-uniform) panels of Fig.^. One characteristic is that density dependence is weak in the 
non- uniform system. This is brought about by a density-dependent increase of Cy„. The reason why non- 
uniformness increases it can be inferred from an expression obtained by a weak-coupling approximation 
for the uniform system 0, 

Cvn ( Ao\^^^ ( Aq 



T 



oc 



T 



exp 



T 



(10) 



with Aq being the uniform gap. This is a decreasing function of ^o/T at low temperatures, Aq/T > 
5/2, and consequently leads to Cy,i (normal) > Cv„ (super). Since non-uniformness produces a region 
where the gap is small, it results in an increase of Cy„. This effect is conspicuous in particular at lower 
temperatures because a small variation in Aq leads to a large one in Ao/T. Conversely, difference in Aq 
is not important at high temperatures. 
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Figure 2: Temperature dependence of the partial specific heat of uniform neutrons (left), non-uniform 
neutrons (center) and electrons (right). 

Since the contribution from electrons is of the same order (the right panel of Fig.^ as that from the 
non-uniform neutrons, the enhancement of Cytotj 

Cytot (nu) _ Cy„ (nu) + Cye ^^^-j 
Cytot (u) C'v„(u) + Cve 

where "nu" and "u" stand for non-uniform and uniform, respectively, is factor 3 - 4 as shown in Fig.^. 
Comparing this to Fig. 3 in Ref. highlights the difference at low temperatures. The consideration above 
indicates that the difference originates from the region where pairing gap is small — the interior region of 
the nucleus at the center of the WS cell. Actually, although the Gogny force adopted in Ref. [Q is known 
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to provide us with pairing properties similar to those given by bare interactions at low densities 
there is a big difference at high densities — superfluid phase survives up to much above the saturation 
density (see Fig. 5. 4 of Ref.^, for example). In contrast, our cutoff parameter in the RMF description 
was determined so as to reproduce the pairing properties given by the Bonn-B potential and accordingly 
pairing gap closes at around the saturation density. This makes the gap in the interior region vanish. In 
Ref. the non-relativistic version of the Bonn-A potential was adopted, and the results for p/ po = 0.058 
and 0.176 in the spherical phase were also presented. Our results resemble theirs. Here one comment is in 
order on the double-counting problem of short-range correlations which arises when effective interactions 
are adopted in the gap equation. In principle this problem exists as pointed out in Ref.Q, but an analysis 
in Ref.H indicates that it is not serious practically. 




T(MeV) 

Figure 3: Temperature dependence of the ratio of the total specific heat. Legend of the curves is the 
same as in Fig.^. 

Since the essential difference between Ref.||l| and ours is the gap in the interior region as mentioned 
above, we examine a schematic interpolation between the two models; we vary the gap in this region by 
hand as shown in the left panel of Fig.||for density (b). The resulting specific heats are presented in the 
center panel. Clearly Cvn approaches to that in Ref.[|l[ (see Fig. 2 there) as the gap in the interior region 
is increased, and to that of the uniform matter by a further increase. Finally the total specific heats are 
shown in the right panel. This confirms that the vanishing gap in the interior region is the very cause of 
the enhancement of Cvtot at low temperatures. 
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Figure 4: Schematic variation of the gap profile (left), resulting temperature dependence of the neutron 
specific heat (center) and of the ratio of the total specific heat (right), at density (b). 

Beyond this first-step calculation, two kinds of refinements are in order; one is to take into account 
non-spherical shapes of lattice nuclei in the high-density part of the inner crust, the other is the so-called 
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proximity effect beyond the local-density approximation, that is, spread of the pair wave function. The 
former was examined in Refs.|^, ^ and shown to strengthen non-uniformness, whereas the latter was 
shown to weaken it . 
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